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Abstract In this paper we consider a random motion of magnetic bright points (MBP) associated 
with magnetic fields at the solar photosphere. The MBP transport in the short time range [0 — 20 
minutes] has a subdiffusive character as the magnetic flux tends to accumulate at sinks of the flow 
field. Such a behavior can be rigorously described in the framework of a continuous time random walk 
leading to the fractional Fokker-Planck dynamics. This formalism, applied for the analysis of the solar 
, subdiffusion of magnetic fields, generalizes the Leighton's model. 
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d , 1 Introduction 

Much of the solar surface phenomena is caused by photospheric convective motions of magnetic flux ele- 
ments. Due to a complicated character of the solar convection, the heliospheric magnetic field contains 
■^j- ' a random component. Leighton (1964) suggested to consider the migration of magnetic regions on the 

£C) , Sun as a simple Brownian random walk on the solar surface. Until the precision of experiments was not 

high, this model was good enough as an approximation to the diffusion of magnetic elements. However, 
the recent experimental studies have shown that the model is too simple to envelop the phenomenon in 
full. The MDI magnetogram data (Hagenaar et al., 1999) from the SOHO spacecraft have noticed that 
the diffusion coefficients in tracking magnetic elements vary in time. Therefore, Cadavid et al. (1999) 
gave a refinement of the Leighton's model. In their examination the walkers (magnetic elements) stick 
before the next jump, i. e. there are traps at stagnation points slowing down the walkers' motion. The 
statistical analysis of the MBP data approves such a subduffisive transport in the intergranular lanes 
within time interval to about 20 minutes. Nevertheless, this point of view produces an impression of 
an unfinished work. The Leighton's model clearly leads to a diffusion equation. It is hence reasonable 



X 

to ask whether it is possible to find a similar equation concerning the subdiffusive case. This question 
remained without any answer in the work of Cadavid et al. (1999). Perhaps, therefore Giacalone 
and Jokipii (2004) stressed that they do not see evidence of anomalous diffusion in random walks 
of magnetic footpoints. However, one may observe that the analysis of Giacalone and Jokipii (2004) 
is related to supergranular scales, with times being the order of days. Their casual comment about 
the results of Cadavid et al. (1999) is somewhat misleading since they did not investigate the same 
phenomena. After all, the analysis of Cadavid et al. (1999) supports the normal diffusion for larger time 
scales. Following Simon et al. (1995), the migration of the magnetic flux from a point to a point tends 
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to accumulate at sinks of the flow field. The sinks displace randomly. Zimbardo, Veltri and Pommois 
(2000) have studied magnetic field line transport in 3D magnetic turbulence with anisotropy in the 
parallel and perpendicular directions to the magnetic field. The transport regime of magnetic field lines 
depends on parameters like the magnetic fluctuation level 5B/B , the correlation lengths of magnetic 
turbulence l x ,ly,lz an d the dimensionality of turbulence. The numerical study (Pommois, Veltri and 
Zimbardo, 1999) shows that the transport can be anomalous (subdiffusive or superdiffusive) . The 
various transport regimes are conveniently classified in terms of the Kubo number (Pommois, Veltri 
and Zimbardo, 2001). For magnetic turbulence this number is defined as R = (5B/Bq)(Iu/1±) , where 
in = l z is the correlation length parallel to the average field Bo = Boe z , l± = l x = l y the correlation 
length perpendicular to Bo. For R < 0.2 there are anomalous non-Gaussian transport regimes, whereas 
for 0.2 < R < 1 there is an approximately quasilinear Gaussian diffusive regime. From this point of 
view the analysis of macroscopic diffusion equations for the motion of magnetic field lines represents a 
great interest in physics and astrophysics. 

Recent progress (Metzler and Klafter, 2000; Zaslavsky, 2002; Stanislavsky, 2004; Meerschaert and 
Schefflcr, 2004; Magdziarz et al., 2007; Magdziarz and Weron 2007) in understanding anomalous dif- 
fusion allows one to represent a subdiffusive transport of magnetic fragments on the solar surface in a 
more comprehensive form. The purpose of this paper is to perform this work. A simple mathematical 
review on the random processes of anomalous diffusion is given in Section 2. This view will be especially 
useful for readers far from the probabilistic theory. In Section 3 we derive the fractional Fokker-Planck 
equation (FPE) in spherical coordinates. This particular equation corresponds to a rotational subdif- 
fusion being of interest to the solar physics of magnetic fields. In next section we discuss the solution 
describing the migration of the MBP on the Sun. Section 5 contains conclusions. 



2 Diffusive processes as a continuous limit of CTRW 

The notion of continuous time random walks (CTRW) has been introduced in physics by Montroll 
and Weiss (1965). They generalized a simple random walk which is based on the assumption that step 
changes (jumps) are made through equal time intervals. In contrast, the CTRW concerns random walks 
with a random waiting time among subsequent random jumps. The generalization has become very 
popular for many physical applications of anomalous diffusion, e. g. for transport in disordered media, 
superslow relaxation, etc. However, this model is useless until random values (jumps and waiting times) 
arc undefined by a probabilistic description (for example, by their probability densities or characteristic 
functions). 

In order to explain the idea of a CTRW let us consider the simplest one-dimensional case. The 
position of a walker after k random jumps is given by 

k 

R(k)=Y, R *> A(0)=0, (1) 

i=l 

where {Ri} are random variables representing the length and the direction of jumps. The corresponding 
time interval reads 

k 

T(fc)=]TT 4 , T(0)=0, (2) 

i=l 

where {T^} represent the random time intervals (waiting time) between successive jumps of a walker. 
The number N t of jumps performed by the walker till time t is a counting process given by the following 
relation 

N t = max{fc : T{k) < t} , (3) 

meaning that the random number Nt of jumps occurred up to time t is equal to the largest index k for 
which the sum T\ +T 2 + . . . Tfc = T(k) of k random time intervals does not exceed time t. Consequently, 
the total distance reached by the walker up to time t determines the stochastic process 



R(N t )=^2Ri, 

i=l 



(4) 
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known as the CTRW. The analysis of probabilistic properties of such a random sum as given in 
(|U) is a core of limit theorems in probability theory. Limit theorem, under a certain necessary and 
sufficient mathematical conditions, yields the continuous limit of the sums (or any other operation 
on the sequences of random variables , following Feller, 1971). In particular, taking into account the 
expected values (Ri) < oo and (Tj) — > oo, we get the continuous limit X(S t ) of the CTRW process 
(H}. Let us note that the infinite expected value of the waiting time Tj is connected with a long-tail 
property of its probability density function: h(t) ~ t _Q_1 for t — > oo and < a < 1 (for more details, 
see Magdziarz and Weron, 2006). The limiting process, known as the anomalous diffusion process, is 
expressed in terms of the parent process X(t) subordinated by the random time clock St- In terms of 
probability density functions we can write the probability density of the subordinated process X(St) 
as an integral relation 

/•oo 

P{x, t)= I f(x,r) g{t,r)dr , (5) 
Jo 

where f(x,r) and g(t,r) are the probability density functions of X and St, respectively. Here, the 
probability density f{x,r) represents the probability of finding the parent process X(t) at x on the 
operational time r, whereas g(t, r) describes the probability for the operational time r to coincide with 
the real time t. The physical interpretation of St is that this process accounts for the amount of time, 
when a walker does not participate in motion (Baeumer et al., 2005). If the walker randomly moves all 
time, the operational time coincides with the physical one, and g(t, r) is simply the Dirac <5-function. 

It is very important that the functional form of the probability density of the random variable St 
can be calculated explicitly. The procedure can be given in few steps. According to Bingham (1971), 
the Laplace transform of the probability density g(t, r) with respect to r equals 



g(t,v) = / e-^g(t,r)dr = {e'^) = E a (~vt a ) 



o 



where E a (—vt a ) is the Mittag-Leffler function. The Mittag-Lefiler function has a simple Laplace image 
with respect to t, namely 

e- ut g(t,v)dt = u a - 1 /(u a + v). 



n 



The latter expression is easily inverted by Laplace with respect to v as an exponential function. As a 
consequence, the Laplace image of g(t, r) with respect to t is 



poo 

/ e- ut g(t,T)dt = u Q - 1 exp{-u Q r}. 



By taking now the inverse Laplace transform of g(u, r) with respect to u, we obtain the probability 
density of the process St in the form 

g(t,r) =r a F a (r/t a ), (6) 
where the function F a (z) can be written as a Taylor series 



i — ' nil 1 — a — not) 



n=0 

Thus, according to Eq. (O, the probability density p(x,t) reads 



3C 



p(x,t) = / F a (z)f(x,t a z)dz (7) 
Jo 

and its Laplace image takes the form 

/•OO 

p(x 7 u)= p(x,t) exp{-ut}di = u a_1 /(x,it Q ), (8) 



o 

where 

/•OO 



/•OO 

/(x,m q )= / /(x,r) exp{-M Q r}dr. 
Jo 
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The representation ([8]) will be used in the next section for derivation of the macroscopic equation of 
anomalous diffusion. 

The above-mentioned mathematical techniques, appropriate for translational motion of a walker, 
can be also applied for a rotational random walk. As for the motion of MBP, associated with magnetic 
fields at the photosphere, subdiffusion on a sphere is of interest to us. This is nothing else but a 
rotational subdiffusion. Then the space jumps will be given by a polar angle 9 and a longitude <f) in 
the spherical coordinates. It is important to observe that the density f(x, r) may be governed by a 
Fokkcr-Plank equation (FPE) with a time-independent potential, discussed in the next section. 



3 Equation of rotational subdiffusion 

To derive an equation of rotational subdiffusion, we start with one-dimensional case. Let L(x) be a 
time-independent Fokker-Planck operator well-known in the classical statistical physics. Assume that 
the probability density f(x, t) describes an ordinary Brownian motion with respect to the operational 
time r. It will satisfy the FPE with the operational time: 

df(x,r)/dr = L(x)f(x,r). 

Applying operator L(x) to the Laplace image ([5]), we find the following expression 

L(x) p(x,u) = u a p(x,u) — q(x) u a ~ , 

where q(x) is the initial condition. The inverse Laplace transform of the latter gives the fractional 
integral form of the FPE 



1 



p(x,t) = q(x) + I dr(t-T) a - 1 L(x)p(x, T ). (9) 
A") Jo 

The kernel of this integral is a power function resulting in long-term memory effects for the process 
X(St). This memory is a direct consequence of subordination of the space variable X(t). For a = 1, 
equation reduces to the ordinary FPE without any memory effects. It should be mentioned that 
the integral form of equation ^ corresponds to the equivalent differential form 

d a P (x,t) g(x)t- a f 

where d a /dt a denotes the Liouville-Riemann fractional differential operator of order a (Samko et al., 
1993). The Liouville-Riemann fractional differential operator of integer order is simply the ordinary 
derivative. 

As a consequence, the equation of rotational subdiffusion in the simplest case of motion on a circle 

d a w a {8,t) w a {e,o)t- a Dg d ( . a dw a (e,ty 



8t a r(l-a) sin 9 89 88 J' ^ 

where Dg is the diffusion coefficient and W a (9,0) the initial condition. Here, the probability density 
W(9,t) depends on the polar angle 9 only. The solution of Eq. fTU)) may be expressed in terms of the 
following integral transformation 



W a (9,t)= F a (z)W 1 (9,t a z)dz. 
Jo 

This implies that the averaged value of cos 9 equals 

(cos6»)= [ cos 6 W a (6,t) sin6 d6 = E a (-2D e t a ) . 
Jo 

For large t all the directions of the walker motion become equiprobable. However, in contrast to the 
normal rotational diffusion this equilibrium is reached slowly because of a power asymptotics of the 
Mittag-LefHer function. 
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4 Transverse subdiffusion of magnetic footpoints at the Sun 

We now write down the equation of subdiffusion associated with the observed motion of magnetic 
footpoints embedded in the transverse flows on the solar surface. The distribution g a (0, cj> ,t) of mag- 
netic field footpoints on the solar spherical surface of radius Rq obeys the two-dimensional equation 
in spherical coordinates 



dt a r(l-a) RQsm9d9\ R Q 89 J R Q sin 9 d<j> \ R Q sin 9 d(j> 

where k is the subdiffussion coefficient with the following dimension [k] = 1 /time a . The solution of 
Eq. can be obtained via a separation ansatz in terms of the spherical harmonics Y mn (9, </>), namely 



oo m—n , >. 

g a {6^,t)=Y J a mn Y mn {9,^)E Q (~-^7i(n + l)t a j. (12) 

n—d m,——n. V © / 



If one takes an impulse release of footpoints at the location (9q, </>o)> the initial condition is determined 
through the Dirac (^-function 

g a (9, 0,0) = -A^ 5{<P - <f> a ) 5(9 - O ). 

Using the completeness relationship of spherical harmonics states, one obtains the amplitude coeffi- 
cients a mn expressed by the complex conjugate of the spherical harmonics Y* ln (9o, </>o)- It is interesting 
to consider an asymptotic behavior of the solution (|12[) as t — >• oo. The distribution approaches asymp- 
totically the limit 

lim g a (9,<p,t) = — , 

t— >oo 47T 

which simply corresponds to the homogeneous distribution on the sphere in equilibrium. After inte- 
grating expression (|T2j) from = to 2ir, we get the solution of Eq. (fTTj) which will depend only on 
the polar angle 9 and time t: 

°° 2n + 1 / k \ 

g a (9,t) = 22^—P n (cos6 )P n (cos9)E a l~— n (n + l)t a \ , (13) 

n=0 V O / 

where P n (y) is the Legendre polynomial of degree n. All the odd moments of the distribution are equal 
to zero, but the even moments are not. So the first even moment takes the form 

<sin^) = -(l-^(-^ 



Since the subdiffusion motion of magnetic footpoints is characterized by a small value k/Rq (about 
0.001 in order of magnitude), we arrive at the relation 



4k 

r{iTa)Rl 



It is obvious that for a = 1 this description leads to the Leighton's model. 

In the framework of the fractional FPE approach to the transport of magnetic fields on the Sun we 
are able to recover the experimental evidence which demonstrates the change of diffusion properties 
during the life time of the MBP, from subdiffusion to normal one. From the theoretical point of view 
the difference between subdiffusion and ordinary diffusion is related with an evolution of the value 
of the parameter a, from a < 1 to a = 1. This change corresponds to different properties of the 
distribution of the interjump time intervals of the MBPs. If the expected value of the waiting time is 
infinite (0 < a < 1), then bright points have to be trapped. This effect is connected with the long- 
tailed properties of the waiting-time distribution, the necessary condition to obtain the subdiffusion 
from the CTRW scheme. For a=l there are no traps, since the expected waiting-time value is finite 
or time variable is deterministic. Due to the traps the diffusion of the MBPs has a mixed character of 
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random stops and motion, whereas in the case of normal diffusion the motion of the MBPs continues 
all time. From the astrophysical point of view the appearance of subdiffusion in the short-time range 
means that the life time of the traps is shorter than the life line of the MBP. Following the magnetic 
turbulence studies on laboratory and astrophysical plasmas, the anomalous diffusion of magnetic field 
lines may be associated with the existence of closed magnetic surfaces. Probably, the surfaces serve as 
traps for the MBP. The slope of the experimental data variance in time is a criterion for revealing the 
peculiarity. The comparison of the subdiffusive model with experimental data may be carried out in 
the same manner as this is the case for the work of Cadavid et al. (1999). 

5 Conclusions 

We have presented an approach to anomalous diffusion which follows from an intuitive concept of 
sticking times suggested in the paper of Cadavid ct al. (1999). We have derived the fractional FPE in 
spherical coordinates appropriate for describing the subdiffusive migration of the MBP. As a special 
case, the normal diffusion of the MBPs on the Sun can be obtained. Therefore, the Leighton's model 
has been generalized. It should be noticed that the consideration of two different diffusive regimes in the 
MBP motion, as resulting from a simple sum of an ordinary first temporal derivative and a fractional 
temporal one in the FPE, will not lead to the expected phenomenon. In this case for short-time region 
such a temporal operator would lead to a normal diffusion, whereas for the large-time scale it does a 
subdiffusion (see, as an example, the work of Schumer et al., 2003). 
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